Tunable thermopower in a graphene-based topological insulator 
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Following the recent proposal by Weeks et al., which suggested that indium (or thallium) adatoms 
deposited on the surface of graphene should turn the latter into a quantum spin Hall (QSH) insulator 
characterized by a sizeable gap, we perform a systematic study of the transport properties of this 
system as a function of the density of randomly distributed adatoms. While the samples are, 
by construction, very disordered, we find that they exhibit an extremely stable QSH phase with no 
signature of the spatial inhomogeneities of the adatom configuration. We find that a simple rescaling 
of the spin-orbit coupling parameter allows us to account for the behaviour of the inhomogeneous 
system using a homogeneous model. This robustness opens the route to a much easier experimental 
realization of this topological insulator. We additionally find this material to be a very promising 
candidate for thermopower generation with a target temperature tunable from 1 to 80 K and an 
efiiciency ZT 1. 



PACS numbers: 

I. INTRODUCTION 

The unique states which exist at the boundary of topo- 
logical phases are of fundamental interest, but also of po- 
tential interest for applications due to their robustness. 
Ideally, these edge states determine entirely the transport 
properties of the sample. Moreover, the topological na- 
ture of the bulk phase implies a robustness of these edge 
states with respect to e.g. sample geometry, impurities 
or other external perturbations. This topological robust- 
ness constitutes a fascinating motivation for designing 
electronic or spintronic devices based on these phases, as 
opposed to more fragile materials like e.g. graphene [I]. 
Of particular interest are the recently discovered topo- 
logical insulators which are the focus of much attention: 
their edge states possess unique spin transport proper- 
ties, due to their fixed spin helicity imposed by the in- 
herent strong spin-orbit (SO) interaction. The surface 
states of the three-dimensional realization of these topo- 
logical insulators resemble the low-energy Dirac fermions 
of graphene but here with a single species (as opposed 
to four for graphene), which implies the expected topo- 
logical robustness • Similarly, the edge states of the 
two-dimensional QSH topological phase are helical, con- 
sisting of a Kramers pair of states carrying two opposite 
spin currents [5j 6j. 

Unfortunately, up to now probing the edge states of 
these new phases through transport experiments has 
proven a major challenge. While the surface states 
of candidate materials for the three-dimensional phase 
have been identified through ARPES and STM experi- 
ments [7HH]) transport experiments have been hampered 
by a large contribution of spurious bulk states which 
obscures any clear signature of the surface states [10]. 
On the other hand, the two-dimensional QSH phase has 
only been discovered experimentally in epitaxially grown 



HgTe/CdTe quantum wells following the initial theoret- 
ical proposal jm I12j . In this material, an unambigu- 
ous signature of the helical edge states was observed 
through fractional four-probe conductance measurements 
(13 . However, only one experimental group has mastered 
the fabrication of HgTe quantum wells so far. 

In this article, we focus on a recent proposal to induce 
the two-dimensional QSH phase on a much more read- 
ily available material: graphene. While the intrinsic SO 
interaction of graphene is too weak to open an observ- 
able gap following the initial Kane and Mele proposal 
[5l|6], it is possible to increase it by extrinsic means, as 
was demonstrated in a recent experiment [14! where the 
generation of very large pseudo-magnetic fields by apply- 
ing mechanical strain was observed. A more recent pro- 
posal suggests that the topological insulating behavior 
can be induced by suitable heavy adatoms deposited on 
the surface of the graphene sheet: uniform coverage with 
indium (or thallium, which is unfortunately poisonous) 
was predicted to induce a QSH phase ^ with an insu- 
lating gap of the order of 100 K for a few percent cov- 
erage by indium. Given the technological expertise that 
exists for graphene, the possibility of a graphene-based 
two-dimensional topological insulator opens formidable 
perspectives. Since the breakthrough scotch-tape exper- 
iment carried by the Manchester group [TB], the field of 
graphene science has grown to impressive proportions, 
largely due to the numerous promises this revolution- 
ary material holds for many applications ranging from 
(opto-)electronics to chemistry [T71IIS]- Graphene sam- 
ples of good quality are now routinely available around 
the world. 

In the present work, we explore the condition and do- 
main of existence of the QSH phase in the realistic sit- 
uation where adatoms are deposited at random on the 
surface of graphene. The signature of the QSH phase 
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is found for an extremely wide range of parameters, sig- 
naling its robustness to inhomogeneous SO interaction. 
More importantly, we show that the transport properties 
of this phase show no sign of this inhomogeneity, and 
can be entirely described by an effective homogeneous 
theory which is valid for any adatom configuration (i.e., 
before and after averaging over adatom configurations). 
We find that the resulting QSH phase shows remark- 
able tunable thermopower characteristics, in a range of 
temperatures where usual (narrow-band semi-conductor) 
materials lead to rather poor yields. 

II. MODEL 

We consider a tight-binding model describing the low- 
energy spectrum of graphene in the presence of adatoms 

{i,j),a y {{i,o))£-y,a,l3 

. (1) 

Indices i,j label the lattice sites, a, /3 the spin quan- 
tum numbers, () stands for nearest neighbors, (()) for 
second nearest neighbors. The nearest-neighbor hopping 
t 3eV will set the energy scale: hereafter all energies 
are expressed in units of t while sizes are expressed in 
units of y/ia (a is the carbon-carbon distance). The sec- 
ond term describes the dominant SO coupling induced 
by an indium (or thallium) adatom, with the usual 
Pauli matrix in spin space. The adatoms are distributed 
at random on a fraction riad of the graphene hexagonal 
plaquettes T (see the lower inset of Fig. [l] for a cartoon of 
a plaquette). SO coupling induces chiral second nearest 
neighbor hoppings with Vij = 1 when moving counter- 
clockwise around a plaquette, and —1 otherwise. The 
intensity of the induced SO interaction is parametrized 
by Ago with Ago ~ 0.0067 for indium adatoms (and 
Aso « 0.02 for thallium). 

The original Kane-Mele Hamiltonian [51 |B] for pure 
graphene is recovered by considering SO coupling on all 
plaquettes (rtad = !)• In this case, the uniform SO in- 
teraction opens a gap 2Aso — 6-\/3Aso in the vicinity of 
the Dirac points if, K' of the Brillouin zone. The cor- 
responding QSH phase possesses a Z2 topological sym- 
metry, which manifests itself through the appearance of 
a pair of edge states at the boundary of a sample. The 
signature of this phase was found through ab-initio cal- 
culations in Ref. |15j . 

III. PARAMETRIC STUDY OF THE QSH 
PHASE 

In order to monitor the existence of the QSH phase, we 
focus on the signature of the associated edge modes on 
the conductance matrix of a multi-terminal sample. The 
latter is given by the multi-terminal Landauer-Biittiker 



formula, which expresses the current la flowing in the 
electrode a (see upper inset of Fig. [I] for a sketch of the 
sample) in presence of a potential Vp as 

Ic.^^^Y.^apVp, (2) 

where Tap is the transmission coefficient between two 
electrodes. It can be computed numerically from the 
knowledge of the retarded Green function G via the for- 
mula 

Tap = Tr[r„Gr^G^] , (3) 

with Ta = Im(I]c<), where Eq. is the self-energy of lead 
a, and with the Green function given by the expression 

G{E) = {E~H~Y.^a)-' (4) 

(with the Hamiltonian H in its first quantization form). 
In the QSH phase the transmission between successive 
probes reads exactly Ta^a+i = 2a,Q-i — Ij while all 
other transmission coefficients Ta^p {a ^ /?) vanish. In 
the following, we will use this unique characterization 
of the QSH phase in the 4-terminal cross geometry of 
Fig. [1] The 4-terminal geometry always provides an un- 
ambiguous characterization of this phase, as opposed to 
the two-terminal conductance where g2T = 2e^//i can be 
observed at the Dirac point independently of the presence 
of SO interaction. 

We also pay special attention to avoid the appearance 
of spurious effects in nano-ribbon geometries. Structural 
boundary conditions of graphene are of two types: arm- 
chair and zigzag. The latter features a zero-energy edge 
state, signatures of which have been shown to obscure the 
appearance of the QSH edge state |19J. Therefore we will 
henceforth consider armchair-terminated graphene sys- 
tems in order to avoid confusion between different types 
of edge states. Armchair ribbons are themselves divided 
into two families [20]: they can be either metallic when 
their width W equals 2 (modulo 3), or semi-conducting 
for other values of the width. In the following, we focus 
on metallic armchair ribbons to avoid the competition 
between the SO-interaction-induced gap and the finite- 
width-induced gap. We use a cross-like geometry, with a 
fixed aspect ratio (see inset of Fig. [T]) and smallest width 
W at contacts 1 and 3: the shape of the sample has been 
chosen such that most of the current is directly trans- 
mitted from to 2 in the absence of SO interaction. The 
numerical calculations were performed using the KNIT 
package which implements a generalization of the recur- 
sive Green function algorithm to multi-terminal systems 

Fig. [1] shows the "longitudinal" Tjo and "Hall" Tio, 
Tso transmission coefficients as a function of the strength 
of the SO interaction and for various concentrations rtad 
of adatoms. One observes, as expected, that upon in- 
creasing the strength of Ago for a fixed concentration riad 



3 




n X 

ad so 

Figure 1: (Color online) : Upper inset: schematic of our setup; 
a 4-terminal graphene cross with armchair edges and width 
of the small arm W = 22. The different colors correspond to 
an actual calculation of the current density for spin-up elec- 
trons upon injection from contact 0. The existence of an edge 
state is manifest. Main figure; Scaling of the transmissions 
Tjo from contact for various couplings Aso and adatom den- 
sities riad, plotted as a function of the effective SO coupling 
strength Ago = Aso^iad- Dashed (black) lines correspond to 
"longitudinal" transmission T20 and dotted lines correspond 
to "Hall" transmissions Tio (red) and T30 (green). The two 
sets of curves correspond to an energy E — (open symbols), 
and E — 0.05 (filled symbols). Different symbol shapes corre- 
spond to different values of SO coupling, Ago = 0.02, 0.05, 0.1 
and 0.15 (because they all collapse on the same curve, various 
symbols may seem indistinguishable from each other). Note 
that no averaging over adatom configurations has been per- 
formed here. Lower inset: cartoon of a graphene plaquette. 
Full lines stand for direct hopping elements while dashed lines 
correspond to the SO-induced hopping elements. 



of adatoms, one enters the QSH phase: the Hall coeffi- 
cients tend to unity while the longitudinal one vanishes. 
The upper inset shows the actual up spin current den- 
sity inside the sample in the QSH regime: we recover the 
expected edge state characteristics of the QSH phase, 
including an exponential decay of the current as a func- 
tion of the distance to the edge but also oscillations with 
frequency |K — K'| coming from the valley- mixing arm- 
chair boundary condition |19) . The crucial point shown 
in Fig. [1] is that all the results are rescaled as a function 
of the effective SO interaction strength 

Aso = AsoHad ■ (5) 

In other words, the QSH phase, while originating from 
a very inhomogeneous sample, is perfectly described in 
each sample by an effective homogeneous phase with a 
uniform coverage but a reduced SO strength. In this 
phase, each adatom occupying one of the plaquettes 
spreads its coupling Ago over a distance ^ad = 1 / V^-ad , 




Figure 2; (Color online); Statistics of the transmission co- 
efficients T20 (lower curve, black) and T^q (upper curve, red) 
when averaged over 2500 adatom configurations, as a function 
of the effective SO coupling strength. The shaded regions rep- 
resent one standard deviation from the mean value which is 
indicated by white symbols (circles for T20 and squares for 
T30). The corresponding density of probability for 3 points 
(green, cyan and violet filled circles, from left to right) are 
given in the inset: as the SO coupling increases, the proba- 
bility distribution P{T2o) shifts from a Gaussian (right green 
and middle cyan) to a log-normal (left violet). Fits are in- 
dicated by dashed lines. All data points were generated for 
an energy E — 0, a, fixed adatom density Wad = 0.1, and the 
same system sizes as in Fig. [l] 



resulting in a weaker effective SO interaction strength 
Ago = Ago /Cad ^ uniform effective full coverage 
riad = 1- The surprising occurrence of this mean- field 
description in each disordered sample and down to very 
small concentrations can be attributed to the large local- 
ization length C of the QSH edge states (see below). This 
large value of ^ also explains why rather large values of 
Ago are required to observe perfect "Hall" transmissions 
in the small samples considered in Figs. [T] and [2j In real, 
/iTO-sized samples, the required value of Ag^ will be much 
smaller since deviations to perfect "Hall" effect are con- 
trolled by the ratio W/^ between the width of the small 
arm W and the localization length ^ (see Fig. [4] below). 

We have found that the above effective homogeneous 
description is very robust and applies for (i) a single real- 
ization of the adatom configuration, i.e. does not require 
disorder averaging, (ii) various sizes from wide samples 
deep into the QSH phase down to narrow samples where 
the edge states on both sides of the sample have a finite 
overlap, (iii) different energies, from the Dirac point up to 
the SO-interaction-induced gap. To illustrate this robust- 
ness, we performed some statistics and show inFig.[2]how 
the transmission coefficients vary from one adatom con- 
figuration to another. The "envelope" curves displayed 
in Fig. [2] correspond to the average transmission plus or 
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Figure 3: (Color online): Transmission coefficients T20 (cir- 
cles) and T30 (squares) as a function of the effective SO cou- 
pling strength, for various values of onsite disorder: V = 
(filled black), V = 0.4 (empty red) and V = 0.8 (empty 
dashed magenta). No qualitative change is brought about 
by disorder, although some deviation from the V = curve 
in the crossover region between the metal and the QSH phases 
starts to be visible at strong enough {V — 0.8) disorder. The 
data were averaged over 50 distinct realizations of disorder 
and adatom configuration for each value of As*. The energy 
was kept at i5 = 0, and the system sizes as in Fig. [l] In- 
set: Longitudinal transmission T20 as a function of energy for 
Ko = (empty symbols) and A^f = 0.02 (filled symbols). 
Circles are for V = and triangles for a single disorder con- 
figuration with V = 0.8. Disorder has no effect on T20 in the 
QSH phase {E < Asa -0.1). 



minus one standard deviation (i.e. for a Gaussian dis- 
tribution there is a 68% probability that the outcome 
will fall inside the envelope for a given adatom configu- 
ration). At small SO coupling, the system is essentially 
ballistic with very small (Gaussian) fluctuations. At large 
SO coupling, the system is deep in the QSH phase with 
small (log-normal) fluctuations (the inset of Fig. [2] shows 
the actual probability distributions). We find that, in 
the crossover between these two limits, the fluctuations 
remain remarkably low for all values of the effective SO 
coupling strength. Additionally, we saw no evidence for 
a breakdown of our mean-field description, even down to 
very low adatom concentrations. 

A celebrated property of a topological phase is its ro- 
bustness with respect to the presence of disorder. The 
persistence of this property in the present context can 
easily be checked by adding on-site disorder (on each site) 
to our Hamiltonian in Eq. ([T]), following the standard pre- 
scription 

Hdis = ^ Vicl^^Ci^a , (6) 

where Vi is a disorder strength randomly distributed in 



the interval [—1^/2,^/2]. The transmission coefficients 
as a function of the effective spin-orbit coupling for sev- 
eral values of V are shown in Fig. |3] As expected for 
a topological phase, the presence of disorder provides no 
qualitative (and hardly any quantitative) modification to 
the general picture described earlier. Deviations from the 
"clean" {V = 0) case remain small, unless the strength 
of disorder reaches (extremely strong) values of the or- 
der of the hopping parameter t. Note that, in the main 
panel of Fig.[3j the Fermi energy is very close to the Dirac 
point. It is well established that disorder has a very small 
effect on the metallic phase (at small Ago) close to the 
Dirac point [551 123] ■ This point can be understood from 
a semi-classical consideration: close to the Dirac point, 
the Fermi wave length diverges so that only very long 
range disorder affects the physics. 

Away from the Dirac point, onsite disorder does affect 
the transport properties of graphene, but not of the QSH 
phase. This is best seen in the inset of Fig. [3] where we 
plot the longitudinal transmission T20 as a function of 
the Fermi energy in the presence/absence of spin-orbit 
coupling and disorder (for a typical sample). When spin- 
orbit coupling is present, adding disorder has no effect 
on T20 so long as E < Ago « 0.1, i.e. inside the QSH 
phase. In the absence of spin-orbit coupling (or above 
the spin-orbit gap), however, disorder strongly affects T2Q 
except at small energies close to the Dirac point. In short, 
we find that the QSH phase is very resilient to onsite 
disorder. 

The above considerations lead us to predict that for the 
present system, the physics of the QSH phase is entirely 
described by the Kane-Mele model (i.e. the full coverage 
case, Had — 1) provided one performs the substitution 
Aso ~^ AsofT-ad- In particular, the known analytical ex- 
pressions for the characteristic scales of the Kane-Mele 
model should apply here. We now explicitly check this 
for the gap Ago [6] of the QSH phase and the width f 
[TO] of the edge states, whose expected expressions read 

A,o = 3V3Af , (7) 

^ 2A,o 12Af ■ ^""^ 

Fig. [4] shows the behavior of the longitudinal transmis- 
sion T20 (averaged over several adatom configurations) 
as a function of the width W. Its exponential decrease 
T20 oc e~^/'> allows for an accurate determination of ^: 
as the width increases, the overlap between two opposite 
edge modes and the associated backscattering decrease. 
An extremely good scaling of the data is obtained with 
an expression of the transmission deduced from Eq. (|8| , 
confirming the above effective homogeneous description. 

As a last test, we now turn to a study of the two- 
terminal transmission T(E) of a rectangular graphene 
ribbon as a function of the Fermi energy Ep. The results 
are shown in Fig. [s] (upper plot): for energies inside the 
topological gap {Ep < Ago), the current is entirely car- 
ried by the edge states leading to T = 2. Above the gap 
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Figure 4: (Color online): Averaged longitudinal transmission 
as a function of W/^ for energy E = 0. Filled symbols: fixed 
Aso ~ 0.1 with nad = 0.2 (circles), 0.15 (squares), 0.1 (tri- 
angles), 0.05 (diamonds). Open symbols: idem but the role 
of Aso and fixed iiad are exchanged, stars: Aso = 0.02. Line: 
Y oc exp — W/^ where ^ is given by Eq. ([8|. Upper inset: 1/^ 
(extracted from the data of the main plot) as a function of 
Aso including additional points for different energy E = 0.02, 
different aspect ratio of the sample, and more values of riad 
and Aso (various symbols). Dashed line: Eq. ([8|. Lower inset: 
Aso as a function of Ago for Aso ~ 0.02, 0.05 and 0.1 (various 
symbols). Dashed line: Eq. In all cases, error bars due 
to sample-to-sample fluctuations are smaller than the symbol 



Figure 5: (Color online): Upper plot: Transmission T as a 
function of the energy _B in a simple armchair ribbon of width 
W, for As^ = 0.01 (left black line: riad = 0.2, crosses: Wad = 
0.5) with W = 130, A^^ = 0.02 (middle magenta hue) with 
W — 65.5, and A^o = 0.04 (right violet line and pluses: two 
difi'erent samples) with W = 32.5. Lower right: ZT as a 
function of temperature 6 for Ago = 0.01, 0.02 and 0.04 from 
right to left. The Fermi energy was chosen to be equal to 
1.5Aso. Lower left: peak value ZTinax of ZT{6) as a function 
of Aso for -^so = 0.02 (diamond), Aso = 0.05 (circles) and 
Aso = 0.1 (crosses). 



[Ep > Ago), one leaves the QSH phase and the transmis- 
sion increases quickly, allowing for a precise extraction of 
the gap value. Note that as Ago gets smaller, wider sam- 
ples must be used to keep the number of open channels 
aX Ep — Ago larger than one. In the lower inset of Fig.|4j 
we plot the extracted Ago as a function of Ag^ for various 
values of riad and find, again, a remarkable agreement 
with the effective homogeneous Kane-Mele description. 



IV. THERMOPOWER 

We now turn to an analysis of a very peculiar feature of 
Fig. [5] the slope of the two-terminal transmission T{E) 
above the gap is very steep, and gets steeper as the gap 
decreases. This property has strong implications in terms 
of thermopower generation. To discuss thermopower at 
finite temperature within the Landauer-Biittiker frame- 
work, let us start by introducing 



dE 



dE 



{E^EfT T{E), 



(9) 



where f{E) = l/(e(^-'^F)/('=B^) + l) is the Fermi function 
at temperature 0, from which we can express the conduc- 
tance g = (e^//i)To, the heat conductance gn = i/{hO)T2 



and the Seebeck conductance gs = —e/{hd)Ti. The 
Seebeck coefficient S = —5V/S9 which characterizes the 
voltage SV across a sample produced by a difference of 
temperature 69 between sample edges is simply given 
by S = gs/g- The dimensionless figure of merit which 
measures the efficiency of the thermopower generator is 
known as the ZT = OS'^u / k parameter (we keep the stan- 
dard notation ZT although we call the temperature 6 and 
not T). a (k) is the electrical (thermal) conductivity. At 
low ZT <C 1, a thermocouple has an efficiency equal to 
the fraction ZT/A of the Carnot efficiency 59/0^ while it 
tends to the Carnot value at large ZT ^ 1. Within the 
Landauer-Biittiker approach, we have |24) 



ZT : 



T0T2 - (Ti)2 



(10) 



which simplifies at low temperature into the so-called 
Mott formula 25 



ZT 



3 ^ ^ ^ V T{Ep) 



(11) 



Hence the key towards an efficient thermocouple lies in 
low values of the transmission and simultaneously a steep 
variation of this transmission with the energy. The shape 
of the transmission versus energy curve T(E) presented 
in Fig. [5] possesses all the required properties: the low 
transmission value T = 2 accounts for the single edge 
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state conducting channel (moreover it does not scale with 
the sample width W) while its derivative is typically of 
order T' « l/Ago- This results in very high values of 
ZT « 1 for temperatures « Aso/fcs- 

To quantitatively describe the expected behavior, wc 
have plotted the parameter ZT as a function of the tem- 
perature 9 in the lower right of Fig. [5j We find well- 
defined maxima for 9 sa Aso/(2A:b), while the optimal ef- 
ficiency is reached for energy values Ep « 3Aso/2. Note 
that the present material possesses the unique character- 
istics that the optimum temperature Aso/(2fcB) can be 
simply tuned by changing the concentration of adatoms 
while the Fermi level can be easily switched to the opti- 
mum value by e.g. a simple back gate. The lower left 
plot of Fig. [5] shows the value of the peak of ZT{9) 
as a function of A^^ : we find very high values, up to 
ZT — 0.5, which tend to increase upon decreasing the 
gap Ago. Hence, these graphene-based topological insula- 
tors appear as very good candidates for low-temperature 
thermocouples: for instance a 6% indium coverage is ex- 
pected to give a gap of 80K, hence an optimum working 
temperature of 40 K. At 1% coverage, the optimum tem- 
perature lies around 10 K, a target temperature for e.g. 
radioisotope thermoelectric generators of spacecrafts for 
which a material with such high ZT w 0.5 would consti- 
tute a great improvement. 

Finally, we note that the above estimate of ZT only 
takes into account the electron contribution to the ther- 
mal conductivity. Real values should renormalize by a 
factor k/{k + Kph) where Kph is the phonon contribution 
and should therefore be slightly smaller than predicted 
above. However, adding structural or extrinsic disorder 
to the ribbon or intentionally damaging the ribbon by 
making holes in it can drastically reduce the phonon- 
mediated thermal conductivity, while the topologically 
protected physics discussed in this paper should remain 
largely unaffected. 

V. CONCLUSION 

One of the most important aspects associated with 
the discovery of graphene is the relative easiness (and 



low cost) of sample production. While very few groups 
can produce high mobility two-dimensional electron gases 
in semiconductor heterostructures, graphene physics re- 
quires lighter equipment and is being studied by an in- 
creasing number of groups. The same remark is even 
more applicable to the study of two-dimensional topo- 
logical insulators. Hence, the proposal [TS] discussed in 
this paper for the realistic situation of inhomogeneous 
samples appears all the more promising. 

We have shown that such a QSH phase, while originat- 
ing from very inhomogeneous samples, can be described 
in each sample by the known results on the pure Kane- 
Mele model with an effective SO coupling strength ac- 
counting for the density of adatoms. This description 
should hold in the presence of disorder or charge-density 
fluctuations induced by the underlying substrate, as long 
as the magnitude of these effects remains smaller than 
the gap Ago- We have further demonstrated that this 
new material and its associated QSH phase provide a 
very efficient thermocouple at low temperature. Indeed, 
while important efforts have been made to improve the 
ZT parameter at room temperature (values of ZT > 1 
can now be found) , the existing materials are very ineffi- 
cient at low temperature. For instance, existing thermo- 
couples working around 10 K have ZT < 0.01. Moreover, 
the graphene-based thermocouples discussed here have a 
target temperature which can be tuned from a few hun- 
dred degrees Kelvin (riad = 1) down to OK (nad = 0) 
by simply changing the concentration of adatoms while 
retaining extremely high values of ZT ss 1. 
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